Contributions of common genetic variants to constitutional delay of puberty and idiopathic hypogonadotropic hypogonadism.

CONTEXT: Constitutional delay of puberty (CDP) is highly heritable, but the genetic basis for CDP is largely unknown. Idiopathic hypogonadotropic hypogonadism (IHH) can be caused by rare genetic variants, but in about half of cases, no rare-variant cause is found. OBJECTIVE: To determine whether common genetic variants that influence pubertal timing contribute to CDP and IHH. DESIGN: Case-control study. PARTICIPANTS: 80 individuals with CDP; 301 with normosmic IHH, and 348 with Kallmann syndrome; control genotyping data from unrelated studies. MAIN OUTCOME MEASURES: Polygenic scores (PGS) based on genome-wide association studies for timing of male pubertal hallmarks and age at menarche (AAM). RESULTS: The CDP cohort had higher PGS for male pubertal hallmarks and for AAM compared to controls (for male hallmarks, Cohen's d = 0.85, p = 1 × 10-16; for AAM, d = 0.67, p = 1 × 10-10). The normosmic IHH cohort also had higher PGS for male hallmarks compared to controls, but the difference was smaller (male hallmarks d = 0.20, p = 0.003; AAM d = 0.10, p = 0.055). No differences were seen for the KS cohort compared to controls (male hallmarks d = 0.04, p = 0.45; AAM d = -0.03, p = 0.86). CONCLUSIONS: Common genetic variants that influence pubertal timing in the general population contribute strongly to the genetics of CDP, weakly to normosmic IHH, and potentially not at all to KS. These findings demonstrate that the common-variant genetics of CDP and normosmic IHH are largely but not entirely distinct.

Selection, optimization and validation of ten chronic disease polygenic risk scores for clinical implementation in diverse US populations.

Polygenic risk scores (PRSs) have improved in predictive performance, but several challenges remain to be addressed before PRSs can be implemented in the clinic, including reduced predictive performance of PRSs in diverse populations, and the interpretation and communication of genetic results to both providers and patients. To address these challenges, the National Human Genome Research Institute-funded Electronic Medical Records and Genomics (eMERGE) Network has developed a framework and pipeline for return of a PRS-based genome-informed risk assessment to 25,000 diverse adults and children as part of a clinical study. From an initial list of 23 conditions, ten were selected for implementation based on PRS performance, medical actionability and potential clinical utility, including cardiometabolic diseases and cancer. Standardized metrics were considered in the selection process, with additional consideration given to strength of evidence in African and Hispanic populations. We then developed a pipeline for clinical PRS implementation (score transfer to a clinical laboratory, validation and verification of score performance), and used genetic ancestry to calibrate PRS mean and variance, utilizing genetically diverse data from 13,475 participants of the All of Us Research Program cohort to train and test model parameters. Finally, we created a framework for regulatory compliance and developed a PRS clinical report for return to providers and for inclusion in an additional genome-informed risk assessment. The initial experience from eMERGE can inform the approach needed to implement PRS-based testing in diverse clinical settings.

Inferring compound heterozygosity from large-scale exome sequencing data.

Recessive diseases arise when both copies of a gene are impacted by a damaging genetic variant. When a patient carries two potentially causal variants in a gene, accurate diagnosis requires determining that these variants occur on different copies of the chromosome (that is, are in trans) rather than on the same copy (that is, in cis). However, current approaches for determining phase, beyond parental testing, are limited in clinical settings. Here we developed a strategy for inferring phase for rare variant pairs within genes, leveraging genotypes observed in the Genome Aggregation Database (v2, n = 125,748 exomes). Our approach estimates phase with 96% accuracy, both in trio data and in patients with Mendelian conditions and presumed causal compound heterozygous variants. We provide a public resource of phasing estimates for coding variants and counts per gene of rare variants in trans that can aid interpretation of rare co-occurring variants in the context of recessive disease.

Pediatric Features of Genetic Predisposition to Polycystic Ovary Syndrome.

CONTEXT: Polycystic ovary syndrome (PCOS) has historically been conceptualized as a disorder of the reproductive system in women. However, offspring of women with PCOS begin to show metabolic features of PCOS in childhood, suggestive of childhood manifestations. OBJECTIVE: To identify childhood manifestations of genetic risk for PCOS. METHODS: We calculated a PCOS polygenic risk score (PRS) for 12 350 girls and boys in 4 pediatric cohorts-ALSPAC (UK), COPSAC (Denmark), Project Viva (USA), and The HOLBÆK Study (Denmark). We tested for association of the PRS with PCOS-related phenotypes throughout childhood and with age at pubarche and age at peak height velocity and meta-analyzed effects across cohorts using fixed-effect models. RESULTS: Higher PRS for PCOS was associated with higher body mass index in midchildhood (0.05 kg/m2 increase per 1 SD of PRS, 95% CI 0.03, 0.07, P = 3 × 10-5) and higher risk of obesity in early childhood (OR 1.34, 95% CI 1.13, 1.59, P = .0009); both persisted through late adolescence (P all ≤.03). Higher PCOS PRS was associated with earlier age at pubarche (0.85-month decrease per 1 SD of PRS, 95% CI -1.44, -0.26, P = .005) and younger age at peak height velocity (0.64-month decrease per 1 SD of PRS, 95% CI -0.94, -0.33, P = 4 × 10-5). CONCLUSION: Genetic risk factors for PCOS are associated with alterations in metabolic, growth, and developmental traits in childhood. Thus, PCOS may not simply be a condition that affects women of reproductive age but, rather, a possible manifestation of an underlying condition that affects both sexes starting in early life.

Phenotyping children and adolescents with obesity using behavioral, psychological, and familial data.

OBJECTIVE: This prospective cohort study aimed to empirically derive phenotypes of children and adolescents with overweight and obesity. METHODS: Latent class analyses using Mplus were carried out in the Growing Up Today Study. Information on participants' weight status, disordered eating behaviors, body image and weight concerns, depressive symptoms, and pubertal timing, as well as and maternal weight status, were included in the latent class analyses, which were stratified by sex. Mixed-effects regression was used to examine associations of the obesity phenotypes with adult weight gain, between age 20 and 35 years, independent of weight at beginning of follow-up and duration of follow-up. RESULTS: Among the girls, four obesity phenotypes were identified: 1) "early puberty"; 2) "mothers with obesity"; 3) "high weight concerns"; and 4) "mixed." Only three phenotypes were identified among the boys: 1) "high weight concerns"; 2) "mothers with obesity"; and 3) "mixed." Participants who had overweight or obesity in childhood or adolescence gained more weight in young adulthood than their leaner peers, but the patterns of weight gain in young adulthood varied by phenotype of obesity in childhood and adolescence. CONCLUSIONS: These results support examining risk factors for and treatment outcomes by obesity phenotypes.

Leveraging polygenic enrichments of gene features to predict genes underlying complex traits and diseases.

Genome-wide association studies (GWASs) are a valuable tool for understanding the biology of complex human traits and diseases, but associated variants rarely point directly to causal genes. In the present study, we introduce a new method, polygenic priority score (PoPS), that learns trait-relevant gene features, such as cell-type-specific expression, to prioritize genes at GWAS loci. Using a large evaluation set of genes with fine-mapped coding variants, we show that PoPS and the closest gene individually outperform other gene prioritization methods, but observe the best overall performance by combining PoPS with orthogonal methods. Using this combined approach, we prioritize 10,642 unique gene-trait pairs across 113 complex traits and diseases with high precision, finding not only well-established gene-trait relationships but nominating new genes at unresolved loci, such as LGR4 for estimated glomerular filtration rate and CCR7 for deep vein thrombosis. Overall, we demonstrate that PoPS provides a powerful addition to the gene prioritization toolbox.

Genetics of skeletal proportions in two different populations.

Human height can be divided into sitting height and leg length, reflecting growth of different parts of the skeleton whose relative proportions are captured by the ratio of sitting to total height (as sitting height ratio, SHR). Height is a highly heritable trait, and its genetic basis has been well-studied. However, the genetic determinants of skeletal proportion are much less well-characterized. Expanding substantially on past work, we performed a genome-wide association study (GWAS) of SHR in ∼450,000 individuals with European ancestry and ∼100,000 individuals with East Asian ancestry from the UK and China Kadoorie Biobanks. We identified 565 loci independently associated with SHR, including all genomic regions implicated in prior GWAS in these ancestries. While SHR loci largely overlap height-associated loci (P < 0.001), the fine-mapped SHR signals were often distinct from height. We additionally used fine-mapped signals to identify 36 credible sets with heterogeneous effects across ancestries. Lastly, we used SHR, sitting height, and leg length to identify genetic variation acting on specific body regions rather than on overall human height.

Prevalence of Deleterious Variants in MC3R in Patients With Constitutional Delay of Growth and Puberty.

CONTEXT: The melanocortin 3 receptor (MC3R) has recently emerged as a critical regulator of pubertal timing, linear growth, and the acquisition of lean mass in humans and mice. In population-based studies, heterozygous carriers of deleterious variants in MC3R report a later onset of puberty than noncarriers. However, the frequency of such variants in patients who present with clinical disorders of pubertal development is currently unknown. OBJECTIVE: This work aimed to determine whether deleterious MC3R variants are more frequently found in patients clinically presenting with constitutional delay of growth and puberty (CDGP) or normosmic idiopathic hypogonadotropic hypogonadism (nIHH). METHODS: We examined the sequence of MC3R in 362 adolescents with a clinical diagnosis of CDGP and 657 patients with nIHH, experimentally characterized the signaling properties of all nonsynonymous variants found and compared their frequency to that in 5774 controls from a population-based cohort. Additionally, we established the relative frequency of predicted deleterious variants in individuals with self-reported delayed vs normally timed menarche/voice-breaking in the UK Biobank cohort. RESULTS: MC3R loss-of-function variants were infrequent but overrepresented in patients with CDGP (8/362 [2.2%]; OR = 4.17; P = .001). There was no strong evidence of overrepresentation in patients with nIHH (4/657 [0.6%]; OR = 1.15; P = .779). In 246 328 women from the UK Biobank, predicted deleterious variants were more frequently found in those self-reporting delayed (aged ≥16 years) vs normal age at menarche (OR = 1.66; P = 3.90E-07). CONCLUSION: We have found evidence that functionally damaging variants in MC3R are overrepresented in individuals with CDGP but are not a common cause of this phenotype.

Corrigendum: Genetic heritability as a tool to evaluate the precision of 24-hour recall dietary questionnaire variables in UK Biobank.

[This corrects the article DOI: 10.3389/fgene.2022.1070511.].

Genome-wide CRISPR screening of chondrocyte maturation newly implicates genes in skeletal growth and height-associated GWAS loci.

Alterations in the growth and maturation of chondrocytes can lead to variation in human height, including monogenic disorders of skeletal growth. We aimed to identify genes and pathways relevant to human growth by pairing human height genome-wide association studies (GWASs) with genome-wide knockout (KO) screens of growth-plate chondrocyte proliferation and maturation in vitro. We identified 145 genes that alter chondrocyte proliferation and maturation at early and/or late time points in culture, with 90% of genes validating in secondary screening. These genes are enriched in monogenic growth disorder genes and in KEGG pathways critical for skeletal growth and endochondral ossification. Further, common variants near these genes capture height heritability independent of genes computationally prioritized from GWASs. Our study emphasizes the value of functional studies in biologically relevant tissues as orthogonal datasets to refine likely causal genes from GWASs and implicates new genetic regulators of chondrocyte proliferation and maturation.

Inferring compound heterozygosity from large-scale exome sequencing data.

Recessive diseases arise when both the maternal and the paternal copies of a gene are impacted by a damaging genetic variant in the affected individual. When a patient carries two different potentially causal variants in a gene for a given disorder, accurate diagnosis requires determining that these two variants occur on different copies of the chromosome (i.e., are in trans) rather than on the same copy (i.e. in cis). However, current approaches for determining phase, beyond parental testing, are limited in clinical settings. We developed a strategy for inferring phase for rare variant pairs within genes, leveraging genotypes observed in exome sequencing data from the Genome Aggregation Database (gnomAD v2, n=125,748). When applied to trio data where phase can be determined by transmission, our approach estimates phase with 95.7% accuracy and remains accurate even for very rare variants (allele frequency < 1×10-4). We also correctly phase 95.9% of variant pairs in a set of 293 patients with Mendelian conditions carrying presumed causal compound heterozygous variants. We provide a public resource of phasing estimates from gnomAD, including phasing estimates for coding variants across the genome and counts per gene of rare variants in trans, that can aid interpretation of rare co-occurring variants in the context of recessive disease.

Epigenome-wide meta-analysis identifies DNA methylation biomarkers associated with diabetic kidney disease.

Type 1 diabetes affects over nine million individuals globally, with approximately 40% developing diabetic kidney disease. Emerging evidence suggests that epigenetic alterations, such as DNA methylation, are involved in diabetic kidney disease. Here we assess differences in blood-derived genome-wide DNA methylation associated with diabetic kidney disease in 1304 carefully characterised individuals with type 1 diabetes and known renal status from two cohorts in the United Kingdom-Republic of Ireland and Finland. In the meta-analysis, we identify 32 differentially methylated CpGs in diabetic kidney disease in type 1 diabetes, 18 of which are located within genes differentially expressed in kidneys or correlated with pathological traits in diabetic kidney disease. We show that methylation at 21 of the 32 CpGs predict the development of kidney failure, extending the knowledge and potentially identifying individuals at greater risk for diabetic kidney disease in type 1 diabetes.

Post-translational control of beige fat biogenesis by PRDM16 stabilization.

Compelling evidence shows that brown and beige adipose tissue are protective against metabolic diseases1,2. PR domain-containing 16 (PRDM16) is a dominant activator of the biogenesis of beige adipocytes by forming a complex with transcriptional and epigenetic factors and is therefore an attractive target for improving metabolic health3-8. However, a lack of knowledge surrounding the regulation of PRDM16 protein expression hampered us from selectively targeting this transcriptional pathway. Here we identify CUL2-APPBP2 as the ubiquitin E3 ligase that determines PRDM16 protein stability by catalysing its polyubiquitination. Inhibition of CUL2-APPBP2 sufficiently extended the half-life of PRDM16 protein and promoted beige adipocyte biogenesis. By contrast, elevated CUL2-APPBP2 expression was found in aged adipose tissues and repressed adipocyte thermogenesis by degrading PRDM16 protein. Importantly, extended PRDM16 protein stability by adipocyte-specific deletion of CUL2-APPBP2 counteracted diet-induced obesity, glucose intolerance, insulin resistance and dyslipidaemia in mice. These results offer a cell-autonomous route to selectively activate the PRDM16 pathway in adipose tissues.

Endocrine and behavioural features of Lowe syndrome and their potential molecular mechanisms.

BACKGROUND: Lowe syndrome (LS) is an X linked disease caused by pathogenic variants in the OCRL gene that impacts approximately 1 in 500 000 children. Classic features include congenital cataract, cognitive/behavioural impairment and renal tubulopathy. METHODS: This study is a retrospective review of clinical features reported by family based survey conducted by Lowe Syndrome Association. Frequency of non-ocular clinical feature(s) of LS and their age of onset was summarised. An LS-specific therapy effectiveness scale was used to assess the response to the administered treatment. Expression of OCRL and relevant neuropeptides was measured in postmortem human brain by qPCR. Gene expression in the mouse brain was determined by reanalysis of publicly available bulk and single cell RNA sequencing. RESULTS: A total of 137 individuals (1 female, 89.1% white, median age 14 years (range 0.8-56)) were included in the study. Short stature (height <3rd percentile) was noted in 81% (n=111) individuals, and 15% (n=20) received growth hormone therapy. Undescended testis was reported in 47% (n=64), and median age of onset of puberty was 15 years. Additional features were dental problems (n=77, 56%), bone fractures (n=63, 46%), hypophosphataemia (n=60, 44%), developmental delay and behavioural issues. OCRL is expressed in human and mouse hypothalami, and in hypothalamic cell clusters expressing Ghrh, Sst, Oxt, Pomc and pituitary cells expressing Gh and Prl. CONCLUSIONS: There is a wide spectrum of the clinical phenotype of LS. Some of the features may be partly driven by the loss of function of OCRL in the hypothalamus and the pituitary.

Orgo-Seq integrates single-cell and bulk transcriptomic data to identify cell type specific-driver genes associated with autism spectrum disorder.

Cerebral organoids can be used to gain insights into cell type specific processes perturbed by genetic variants associated with neuropsychiatric disorders. However, robust and scalable phenotyping of organoids remains challenging. Here, we perform RNA sequencing on 71 samples comprising 1,420 cerebral organoids from 25 donors, and describe a framework (Orgo-Seq) to integrate bulk RNA and single-cell RNA sequence data. We apply Orgo-Seq to 16p11.2 deletions and 15q11-13 duplications, two loci associated with autism spectrum disorder, to identify immature neurons and intermediate progenitor cells as critical cell types for 16p11.2 deletions. We further applied Orgo-Seq to identify cell type-specific driver genes. Our work presents a quantitative phenotyping framework to integrate multi-transcriptomic datasets for the identification of cell types and cell type-specific co-expressed driver genes associated with neuropsychiatric disorders.

Genome-wide meta-analysis and omics integration identifies novel genes associated with diabetic kidney disease.

AIMS/HYPOTHESIS: Diabetic kidney disease (DKD) is the leading cause of kidney failure and has a substantial genetic component. Our aim was to identify novel genetic factors and genes contributing to DKD by performing meta-analysis of previous genome-wide association studies (GWAS) on DKD and by integrating the results with renal transcriptomics datasets. METHODS: We performed GWAS meta-analyses using ten phenotypic definitions of DKD, including nearly 27,000 individuals with diabetes. Meta-analysis results were integrated with estimated quantitative trait locus data from human glomerular (N=119) and tubular (N=121) samples to perform transcriptome-wide association study. We also performed gene aggregate tests to jointly test all available common genetic markers within a gene, and combined the results with various kidney omics datasets. RESULTS: The meta-analysis identified a novel intronic variant (rs72831309) in the TENM2 gene associated with a lower risk of the combined chronic kidney disease (eGFR<60 ml/min per 1.73 m2) and DKD (microalbuminuria or worse) phenotype (p=9.8×10-9; although not withstanding correction for multiple testing, p>9.3×10-9). Gene-level analysis identified ten genes associated with DKD (COL20A1, DCLK1, EIF4E, PTPRN-RESP18, GPR158, INIP-SNX30, LSM14A and MFF; p<2.7×10-6). Integration of GWAS with human glomerular and tubular expression data demonstrated higher tubular AKIRIN2 gene expression in individuals with vs without DKD (p=1.1×10-6). The lead SNPs within six loci significantly altered DNA methylation of a nearby CpG site in kidneys (p<1.5×10-11). Expression of lead genes in kidney tubules or glomeruli correlated with relevant pathological phenotypes (e.g. TENM2 expression correlated positively with eGFR [p=1.6×10-8] and negatively with tubulointerstitial fibrosis [p=2.0×10-9], tubular DCLK1 expression correlated positively with fibrosis [p=7.4×10-16], and SNX30 expression correlated positively with eGFR [p=5.8×10-14] and negatively with fibrosis [p<2.0×10-16]). CONCLUSIONS/INTERPRETATION: Altogether, the results point to novel genes contributing to the pathogenesis of DKD. DATA AVAILABILITY: The GWAS meta-analysis results can be accessed via the type 1 and type 2 diabetes (T1D and T2D, respectively) and Common Metabolic Diseases (CMD) Knowledge Portals, and downloaded on their respective download pages ( https://t1d.hugeamp.org/downloads.html ; https://t2d.hugeamp.org/downloads.html ; https://hugeamp.org/downloads.html ).

Genetic heritability as a tool to evaluate the precision of 24-hour recall dietary questionnaire variables in UK Biobank.

A variety of statistical approaches in nutritional epidemiology have been developed to enhance the precision of dietary variables derived from longitudinal questionnaires. Correlation with biomarkers is often used to assess the relative validity of these different approaches, however, validated biomarkers do not always exist and are costly and laborious to collect. We present a novel high-throughput approach which utilizes the modest but importantly non-zero influence of genetic variation on variation in dietary intake to compare different statistical transformations of dietary variables. Specifically, we compare the heritability of crude averages with Empirical Bayes weighted averages for 302 correlated dietary variables from multiple 24-hour recall questionnaires in 177 K individuals in UK Biobank. Overall, the crude averages for frequency of consumption are more heritable than their Empirical Bayes counterparts only when the reliability of that item across questionnaires is high (measured by intra-class correlation), otherwise, the Empirical Bayes approach (for both unreliably measured frequencies and for average quantities independent of reliability) leads to higher heritability estimates. We also find that the more heritable versions of each dietary variable lead to stronger underlying statistical associations with specific genetic loci, many of which have well-known mechanisms, further supporting heritability as an alternative metric for relative validity in nutritional epidemiology and beyond.

Evidence From Men for Ovary-independent Effects of Genetic Risk Factors for Polycystic Ovary Syndrome.

CONTEXT: Polycystic ovary syndrome (PCOS) is characterized by ovulatory dysfunction and hyperandrogenism and can be associated with cardiometabolic dysfunction, but it remains unclear which of these features are inciting causes and which are secondary consequences. OBJECTIVE: To determine whether ovarian function is necessary for genetic risk factors for PCOS to produce nonreproductive phenotypes. DESIGN, SETTING, AND PARTICIPANTS: Cohort of 176 360 men in the UK Biobank and replication cohort of 37 348 men in the Estonian Biobank. MAIN OUTCOME MEASURES: We calculated individual PCOS polygenic risk scores (PRS), tested for association of these PRS with PCOS-related phenotypes using linear and logistic regression and performed mediation analysis. RESULTS: For every 1 SD increase in the PCOS PRS, men had increased odds of obesity (odds ratio [OR]: 1.09; 95% CI, 1.08-1.10; P = 1 × 10-49), type 2 diabetes mellitus (T2DM) (OR: 1.08; 95% CI, 1.05-1.10; P = 3 × 10-12), coronary artery disease (CAD) (OR: 1.03; 95% CI, 1.01-1.04; P = 0.0029), and marked androgenic alopecia (OR: 1.03; 95% CI, 1.02-1.05; P = 3 × 10-5). Body mass index (BMI), hemoglobin A1c, triglycerides, and free androgen index increased as the PRS increased, whereas high-density lipoprotein cholesterol and SHBG decreased (all P < .0001). The association between the PRS and CAD appeared to be completely mediated by BMI, whereas the associations with T2DM and marked androgenic alopecia appeared to be partially mediated by BMI. CONCLUSIONS: Genetic risk factors for PCOS have phenotypic consequences in men, indicating that they can act independently of ovarian function. Thus, PCOS in women may not always be a primary disorder of the ovaries.

Growth Hormone Stimulation Testing Patterns Contribute to Sex Differences in Pediatric Growth Hormone Treatment.

INTRODUCTION: Males are twice as likely as females to receive pediatric growth hormone (GH) treatment in the USA, despite similar distributions of height z (HtZ)-scores in both sexes. Male predominance in evaluation and subspecialty referral for short stature contributes to this observation. This study investigates whether sex differences in GH stimulation testing and subsequent GH prescription further contribute to male predominance in GH treatment. METHODS: Retrospective chart review was conducted of all individuals, aged 2-16 years, evaluated for short stature or poor growth at a single large tertiary referral center between 2012 and 2019. Multiple logistic regression models were constructed to analyze sex differences. RESULTS: Of 10,125 children referred for evaluation, a smaller proportion were female (35%). More males (13.1%) than females (10.6%) underwent GH stimulation testing (p < 0.001) and did so at heights closer to average (median HtZ-score -2.2 [interquartile range, IQR -2.6, -1.8] vs. -2.5 [IQR -3.0, -2.0], respectively; p < 0.001). The proportion of GH prescriptions by sex was similar by stimulated peak GH level. Predictor variables in regression modeling differed by sex: commercial insurance predicted GH stimulation testing and GH prescription for males only, whereas lower HtZ-score predicted GH prescription for females only. CONCLUSIONS: Sex differences in rates of GH stimulation testing but not subsequent GH prescription based on response to GH stimulation testing seem to contribute to male predominance in pediatric GH treatment. That HtZ-score predicted GH prescription in females but not males raises questions about the extent to which sex bias - from children, parents, and/or physicians - as opposed to objective growth data, influence medical decision-making in the evaluation and treatment of short stature.

Genes with specificity for expression in the round cell layer of the growth plate are enriched in genomewide association study (GWAS) of human height.

Human adult height reflects the outcome of childhood skeletal growth. Growth plate (epiphyseal) chondrocytes are key determinants of height. As epiphyseal chondrocytes mature and proliferate, they pass through three developmental stages, which are organized into three distinct layers in the growth plate: (i) resting (round), (ii) proliferative (flat), and (iii) hypertrophic. Recent genomewide association studies (GWASs) of human height identified numerous associated loci, which are enriched for genes expressed in growth plate chondrocytes. However, it remains unclear which specific genes expressed in which layers of the growth plate regulate skeletal growth and human height. To connect the genetics of height and growth plate biology, we analyzed GWAS data through the lens of gene expression in the three dissected layers of murine newborn tibial growth plate. For each gene, we derived a specificity score for each growth plate layer and regressed these scores against gene-level p values from recent height GWAS data. We found that specificity for expression in the round cell layer, which contains chondrocytes early in maturation, is significantly associated with height GWAS p values (p = 8.5 × 10-9 ); this association remains after conditioning on specificity for the other cell layers. The association also remains after conditioning on membership in an "Online Mendelian Inheritance in Man (OMIM) gene set" (genes known to cause monogenic skeletal growth disorders, p < 9.7 × 10-6 ). We replicated the association in RNA-sequencing (RNA-seq) data from maturing chondrocytes sampled at early and late time points during differentiation in vitro: we found that expression early in differentiation is significantly associated with p values from height GWASs (p = 6.1 × 10-10 ) and that this association remains after conditioning on expression at 10 days in culture and on the OMIM gene set (p < 0.006). These findings newly implicate genes highlighted by GWASs of height and specifically expressed in the round cell layer as being potentially important regulators of skeletal biology. © 2021 American Society for Bone and Mineral Research (ASBMR).